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The construction of Dominic Joyce referred to in the title is for self-dual 
Riemannian metrics with two Killing fields [5]. In this note we show that this 
construction (in the form described by Calderbank and Pedersen in [5]) can be 
extended to a family of nonlinear fourth order PDE in two dimensions, essen- 
tially reducing them to linear equations. One equation in this family is the 
affine maximal equation, and we will see that the construction in this case is 
amost the same as one due to Chern and Terng. 

1 The main result 

We consider convex functions u(xi) defined on a domain in R™ and write J = 
det(iijj), where (wy) is the Hessian ( g§7g^. ) ■ Let ip be any smooth, strictly 
convex function on the half-line (0,oo) and consider the functional 



where (u y ) is the matrix inverse of (ity). This a nonlinear fourth order PDE 
for the function u. In the case when 4>(J) = —J a for some a £ (0, 1) these 
equations have been studied by Trudinger and Wang [6j. The case covered by 
Joyce's original construction is when n = 2 and ip(J) = —log J, as we will 
discuss further in Section 3. (Of course we only consider the functional J 7 as a 
motivation for writing down the partial differential equations (2), and the actual 
convergence of the integral (1) is irrelevant.) Let us say that a point (x\, X2) is 
an "ordinary point" if the derivative V J does not vanish there. 
To state our result, let / be a solution of the equation 




(1) 



The corresponding Euler-Lagrange equations 5T = arc 





/'(*) = ivy ( i). 



Thus /' > and / is a 1 — 1 map from (0, oo) to some finite or infinite interval 
I. Let p be the inverse map raised to the power —1/2, so J -1 / 2 = p(r) when 
r = f(J). Now consider the linear, second order PDE for a function £(H,r) 
defined on some domain in R x /; 

dH ■ 1 df ( 3) 



dH 2 p(r) dr \ dr 

Theorem 1 Suppose £1,^2 are two solutions of equation (3) and dct ( j^jj^y^ 
is positive at some point (H ,ro). Define three 1-forms ei,e 2 ,e by 

ei = p ^f dH ~ §*) - - ^ J f dH + M dr) 

e = £id + 6^2- 

Then e, £1,62 are closed 1-forms and e\ A £2 is non-vanishing near (Ho,ro). 
Thus we can find functions u,xi,X2 with du = e,dxi — ei,dx2 = £2 and xi,X2 
give local co-ordinates around (H ,ro)- If we regard u as a function o/(xi,a; 2 ) 
it is a solution of the equation (2), on a suitable domain. All points in this 
domain are ordinary points. Conversely, any solution of (2) with n = 2, in the 
neighbourhood of an ordinary point, is obtained in this way, with solutions £1,^2 
of (3) which are unique up to translation in the H -variable and the addition of 
constants 

In sum, the local study of the nonlinear equation (2), in two dimensions, is 
essentially equivalent to that of the linear equation (3). 



2 The proof 

This is entirely elementary. From now on we always suppose the dimension n is 
2, and we work locally so we will not specify the precise domains of definition of 
the various functions. Given a convex function u we consider the Riemannian 
metric g — ^Uijdxidxj, and in particular the conformal structure which this 
defines. Write & for the Legcndrc transform co-ordinates = J^-. Suppose now 
that we have some other local co-ordinates Ai, A 2 , inducing the same orientation 
as £1,2:2 (i.e. det (ff 1 ) > 0). We can write the metric g = g a bd\ a dXb in 
these co-ordinates. Recall that the co-ordinates are called isothermal if the 
matrix (g a b) is a multiple of the identity matrix at each point, or in other words 

g = V(d\j + d\l), 

for a positive function V(\i,\z). Let ejj denote the alternating tensor £12 = 
—£21 = 1, en = £22 = 0. 
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Observation 1 The co-ordinates Ai, A2 are isothermal if and only if the partial 
derivatives ^p- , are related by the four equations 

where J = J(x(Ai, A2)). 
To see this we write 

9 = ^ Uijdxidxj = ^ dijdxj = ^ ff^^" rfAarfAfc ' 

«J 3 j-a.b a 



So the isothermal condition is 

d\ a dX h 



ST 9 & dx i 
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In matrix not at iou if .1 = ( ff^ , B = (^§\^j > this i s 



or in other words 



B T = VA- 1 . (5) 



Taking determinants we have dct B dct A = V 2 . On the other hand the matrix 
(mj) is so J — detSdet A" 1 . Thus V = VJdetA and (5) is 

B T = VJdebA A' 1 , 
which is the same as (4) , by the formula for the inverse of a 2 x 2 matrix. 



Now consider any positive function P(X 1 , A 2 ) and the pair of second order 
linear PDE 

Observation 2 7/^(Ai,A 2 ) is a solution of (6) then there is a solution x{X\, A2) 
of the first order system 

d\a dX b 

unique up to the addition of a constant, and x(Xi, X2) satisfies (7). 





d 


+ dX 2 








+ dX2~ 
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This is because the consistency condition for the first order system is 

JL ( p dC\ = _d_ (_ p dt_ 

d\ 2 V dX 2 J <9Ai V 9Ai 

which is the equation (6). In a different language, we are saying that the 1-form 
e = P _1 || 2 (£\i — PJj^d\ 2 is closed, so can be written as dx, for a function x. 
Note that, changing P to there is a complete symmetry between x and £ 

so we can also start with a solution of (7) and construct a solution of (6). 

Observation 3 Suppose £i,£ 2 are two solutions of (6), with det (jjx^j > 0- 
Let x 2 be the solution of (7) corresponding to £1 by (8), and x\ be the solution 
corresponding to — £ 2 . Then det ^ff 1 ) > 0, so xi,x 2 give local co-ordinates. 
Write 

e = £idxi + £2^2 • 

Then e is a closed 1-form and so e — du for a function u. If we express u as a 
function of xi,x 2 then J^- = £3. 

This is all straightforward. The conditions (8) imply that 
det(^WW!^>0. 



\d\aj ' \d\ a 

We have de — d^idxi + d£ 2 dx 2 and, writing £ ija = J^, 

de = (a,irfAi+a i 2rfA2)(-6,2rfAi+6,irfA2)+(6,irfAi+6,2rfA2)(^i,2rfAi-a i irfA 2 ) = 0. 



Now return to our functions ip(J), f(J) and the Euler-Lagrange equation 

(2). 

Observation 4 A convex function u satisfies equation (2) if and only if f(J) 
is harmonic with respect to the metric g = ^Uijdxidxj . 

This is true in any dimension. The formula for the derivative of an inverse 
matrix is 

® ij _ ip 
— — U — — 2_^ U U pqkU q j, 

pq 

whereas the formula for the derivative of the determinant is 

dJ _ , \ pq 
g x . ~ J 2-~, u u pv 

1 pq 

These yield the identity 

E ^ = -E »*W = -E « <p J "|f • w 

j 3 pq j p P 
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Our Euler-Lagrange equation (2) is 

dxi 



where 



3 J 



So 



By the definition of the function / we have 



Using (9) we obtain 



dxj 



Thus the equation (10) is the Laplace equation in the metric g: 



With these four observations the main result, Theorem 1, is almost obvious. 
Suppose we start with an ordinary point of a solution u to (2). Then r = 
/(J) is harmonic, by Observation 4, and we can suppose that its derivative 
does not vanish in the region considered. There is then a conjugate harmonic 
function H , which by definition is one such that the local co-ordinates (H, r) are 
isothermal. By Observations 1 and 2 the functions Xi,£j satisfy the equations 
(6), (7) respectively with Ai = H, A2 = r and P = p(r). Since P docs not 
depend on H the equation (6) can be written in the form (3). By Observations 
2 and 3 we can recover the original function u from the two solutions £ij£;2 
of the linear PDE (or, equally well, the two solutions X\,X2)- It is also clear 
that, conversely, starting with any two solutions £1, £2 to the linear equation we 
construct a solution to (2) by this method. 

Note that the conjugate function H, in this situation, can be defined simply 
as the solution of the system 



dH 

dxi 



= eijVj, (11) 



so we can also thing of H as the Hamiltonian generating the area-preserving 
vector field v in the (#1,0:2) plane. We leave this as an exercise for the reader. 
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3 Examples and discussion 



1. Let , £2) be the Legendre transform of u{x\,X2)- The Hessian ( ^ ^ ) 

at the point £i = J^- is the inverse of the Hessian of u at i,. Thus we can 
write 

where J» = det( q^.q^ )■ This implies that the Legendre transform takes 
a solution u of the equation (2) associated with ip to a solution u* of the 
equation associated to the function 

In our construction this just means replacing the function p(r) by p(— r) _1 
and interchanging the roles of the co-ordinates Xi, £j. 

2. If = -log* the equation (2) is J2 a^-al ■ = °- This is the Ra- 
tion defining a zero scalar curvature Kahler metric in "symplectic" co- 
ordinates, see PP, [?]. Explicitly we introduce two further co-ordinates 
Ox, 62 and consider the Riemannian metric, in four dimensions, 

u^dxidxj + u 1 -' dOidQj. 

The well-known fact that, in two complex dimensions, such metrics are 
self-dual gives the link with Joyce's original formulation of his construc- 
tion. Under the Legendre transform we get another description corre- 
sponding to the function ip*(t) = tlogt. This leads to the equations 
describing zero scalar curvature metrics in "complex" co-ordinates. When 
ijj(t) = — logi we get p(r) — r and equation (3) is the familiar equation 
denning axi-symmetric harmonic functions on R 3 . 

3. For any function p(r) there is an obvious solution £ = H to (3). Thus we 
get a special family of solutions to (2) with £1 = H and £2 some other 
solution of (3). There is another special family, which corresponds to this 
under the Legendre transform, when £1 is a function of r only, so xi = H . 
These correspond to second order equations of Monge- Ampere type which 
give special solutions of (2). For example, in the zero scalar curvature 
case above we have the special solutions where £1 = logr. The function 
u satisfies the equation det ( gj^rr ) = e~^/ 2 and the corresponding four 
dimensional Riemannian metric is Ricci-flat. 

4. For the Trudinger-Wang equations, when ip(t) = —t a with < a < 1/2, 
we get p(r) = rV( 1— 2o (up to factor which is irrelevant because it drops 
out of the equation (3)). If a = 1/2 we get p(r) = e r l 2 . 
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5. Consider the graph of a function u(x\, x%) as a surface in R 3 . The Gauss 
curvature is 

det(uij) J 
(1 + |V«| 2 ) 2 (l + |Vu| 2 ) 2 ' 

and the induced area form is 

dA=(l+ \Vu\ 2 ) 1/2 dxidx 2 . 

Thus 

K 1/4 dA = J 1/4 d Xl dx 2 . 

The left hand side is invariant under Euclidean transformations of R 3 
while the right hand side is invariant under unimodular affine transfor- 
mations of R 2 . Since these two groups generate all unimodular affine 
transformations of R 3 we see that this 2-form is an affine invariant of the 
surface. The graphs of the solutions of equation (2) when tp(t) = t 1 / 4 are 
"affine maximal" surfaces in R 3 . According to Chern and Terng [3j these 
surfaces can be described locally as follows. Let F\ 1 F2,Fz be harmonic 
functions of variables Ai, A2 and write F = F_(Xi, A2) for the correspond- 

0^ F 3^ F 

ing vector-valued function. Then the condition -^5= + -7^5= = implies the 
consistency of the first order system 

dZ dF dZ dF 

0A1 0A2 aA2 0A1 

so there is a solution Z_ which, under appropriate non-degeneracy condi- 
tions, parametrises a surface in R 3 . Chern and Terng show that these 
surfaces are precisely the affine maximal surfaces. We want to relate this 
description to ours. Notice first that our equations can be written in a 
similar form. Given a pair of solutions £1 , £2 of (3) we define a vector 
valued function 5 with components £1,^2, 1. Then if Z_ = (x\,x-2,u) our 
description is the first order system 

dZ . . „ d£, dZ . . Q £ 

In the case when ip(t) = t 4 ! 4 we get p(r) — r 2 so the equation (3) is not 
the ordinary Laplace equation, in the variables (H, r) . However given a 
function F(H,r) we consider the first order system 

, n d£ dF , i( 9F 
p{r) dH =r dH p{r) d-r =r ^- R (W) 

If p(r) = r 2 one readily checks that this system is consistent, so there is 
a solution £ (H, r) . The equation (3) for £ is equivalent to the ordinary 
Laplace equation for F. Now take two harmonic functions F\, F% and set 
Fs{H, r) — r. Then a few lines of calculation show that the system (12) is 
identical to the system (13), when in the latter we use the functions £1, £2 
corresponding to F\,Fz by solving (14). 
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